Onset of Wave Drag due to Generation of Capillary-Gravity Waves by a Moving 

Object as a Critical Phenomenon 
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The onset of the wave resistance, via generation of capillary gravity waves, of a small object 
moving with velocity V, is investigated experimentally. Due to the existence of a minimum phase 
velocity V c for surface waves, the problem is similar to the generation of rotons in superfluid helium 
near their minimum. In both cases waves or rotons are produced at V > V c due to Cherenkov 
radiation. We find that the transition to the wave drag state is continuous: in the vicinity of the 
bifurcation the wave resistance force is proportional to yV — V c for various fluids. 

PACS numbers: 47.35. +i, 68.I0.-m 



An object, moving uniformly on a free fluid surface, is 
subjected to a drag force, that can be of different physical 
nature. The most common one is a viscous drag which 
at low Reynolds numbers, Re < 1, is the Stokes drag, 
F v — Snr/dV, that is proportional to an object velocity, 
V, and at higher Re is the drag which originates from ei- 
ther laminar or turbulent wakes |IJ . (Here d is the sphere 
diameter, and rj is the viscosity.) However, there exists 
another drag force, which originates from generation of 
capillary-gravity waves by an uniformly moving object. 
These waves remove momentum from the object into in- 
finity, that produces the wave resistance force, acting on 
the object 

The dispersion relation for the surface waves is jlj 



= gk + ok? I p, 



(1) 



where u> is the circular frequency, k is the wave number, p 
is the fluid density, g is the gravity acceleration, and a is 
the surface tension. According to Eq. (1), a phase veloc- 
ity of the waves c = to/k has a mimimum V c = (^ga / p) 1 ^ 4 
at the capillary wave number k c — \J pg/ 1 Below V c the 
surface waves cannot be emitted, so no wave resistance 
force acts on the object. The stationarity condition of 
the wave pattern in the frame moving with the object, 
leads to [01 



lo = kVcosO, 



(2) 



where 9 is the angle between V and k. 

Eliminating lu from equations (1) and (2), one finds 



cos6>(fc) = c{k)/V. 



(3) 



This equation has evidently no solutions for V < V c and 
describes the opening of the Cherenkov radiation cone 
at V > V c ||. Generation of surface waves is analogous 
to the Cherenkov emission of electromagnetic waves by 
a charge uniformly moving with a superluminal velocity 
in a medium. Due to the linear dispersion relation the 
Cherenkov radiation of the electromagnetic waves occurs 



at every wave number into a cone with an opening 9 de- 
fined by Eq. (3). 

For a given V > V c the system of equations (1) - (2) 
has solutions in the range > 9 < x, where the angle \ 
is denned by cosx = V c /V. This range of 9 corresponds 
to the interval k\ < k < k^ of wave numbers, resulting 
from the intersection of Eqs. (1) with a straight line of 
slope V. 

Therefore, unlike the electromagnetic waves, the 
capillary-gravity waves are radiated within a finite range 
Afc = &2 — It tends to zero as Afc = 4fc cV / e, when 
V approaches V c from above, where e = (V — V c )/V c . 
The presence of the minimum on the dispersion curve of 
the capillary-gravity waves (eq.(l)) results in a gap in the 
spectrum w(fc), similar to the gap in the energy spectrum 
of a superfluid Q . The existence of the minimum phase 
velocity for the dispersion curve (1) is analogous to the 
existence of the Landau critical velocity for the phonon- 
roton energy spectrum in a superfluid helium Q. The 
problem of the drag onset at the Landau critical velocity 
Vo due to the radiation of rotons was recently considered 
theoretically in Ref. || . Pomeau and Rica observed the 
onset of the roton generation at a certain critical veloc- 
ity, in numerical simulations of the generalized nonlinear 
Shrodinger equation, and pointed out but without proof, 
that the roton drag force close to the onset depends on 
the control parameter e = (V — Vq)/Vq as y/e in 3D and 
as log e in 2D. 

The onset of the wave resistance and its behavior in 
the supercritical region for the capillary-gravity surface 
waves was theoretically considered in recent papers |p). 
The main result of the calculations is the prediction of 
a discontinuous transition to a nonzero wave resistance 
state at V = V c , if the object size is much smaller than 
the capillary wavelength A c = 2it/k c ||. 
In this Letter we investigate the dependence of the wave 
drag force, on the object velocity, V, in the vicinity of the 
transition for various fluids and object sizes. The results 
of the experiments clearly demonstrate that the transion 
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to the wave drag state is continuous one in contradiction 
with the theoretical prediction |^. 

The experiments were conducted in a circular rotat- 
ing channel of 12.4 cm outer radius and 3.6 cm width, 
made of plexiglass, with a stainless steel ball, as an ob- 
ject, suspended on a wire, which is glued to an elastic 
fiber. The channel was driven by a stepping motor via 
a belt transmission with ratio 0.138 and with a velocity 
controlled better than 0.02%. Balls of various diame- 
ters, d, (1.57, 2.35, and 3.14 mm) and wires (0.3 and 
0.7 mm dia) half- immersed into a fluid, were used. Each 
of them is much smaller than A c for any fluid used in 
the experiments. A drag force as a function of the ve- 
locity was measured in two ways. The first method was 
to measure the drag force, applied to an object, when 
the object deviates from its equilibrium vertical position 
with increasing driving velocity in accordance with the 
drag force. The second method was to measure the same 
force but with an object, hold in an initial vertical posi- 
tion by using an appropriate feedback loop. The latter 
was actuated via eddy current commercial displacement 
gauge (EMD 1050 Electro Corporation), that measures 
a deviation of brass cylinder(2 mm dia), mounted on the 
same wire as the object, from its initial position. A signal 
from the displacement detector is used via an instrumen- 
tational amplifier and a power amplifier for the feedback 
loop to drive electromagnets which hold the object in its 
initial position (see Fig.l). A current supplied to the 
electromagnets is calibrated as a force. In another setup 
a ball deviation from its initial position was measured op- 
tically by CCD camera with spatial resolution of about 
10 fim. Light refraction visualization was also used to de- 
fine the onset of the surface waves generation, the angle 
of the Cherenkov cone opening, and the wavenumbers as 
a function of velocity. As working fluids water, silicone 
oils DC200 of the kinematic viscosities v = r]/p =10 and 
50 cS and of the same surface tension a —21.2 mN/m, 
water-glycerol mixtures of v =18, 30, and 46 cS and al- 
most the same a =66 mN/m, were used. The depth of 
each fluid was sufficiently large to ensure the validity of 
an infinite depth approximation for the surface waves. 

The raw data on the full drag force, F, as a function 
of the channel velocity for silicone oil(10 cS)-air inter- 
face and d =3.14 mm ball are presented in Fig. 2. The 
data taken with increasing and decreasing values of veloc- 
ities are reproducible and show no hysteresis. Analogous 
data for d =2.35 and d =1.57 mm balls were also ob- 
tained. Since the critical velocity of the transition to the 
wave generation state is rather low, the critical Reynolds 
numbers at the transition for different fluids and ball sizes 
are in the range between 2.2 (for 1.57 mm ball and sil- 
icone oil 50 cS) and about 700 (for 3.14 mm ball and 
water). It means that below the transition the drag con- 
sists of a viscous Stokes drag at Re < 1, which changes 
linearly with V , and the drag, proportional to V 2 , when 
a boundary layer still remains laminar up to Re ~ 10 3 



|l[ . One fits the data below the transition by a second or- 
der polynomial and gets a reasonable value of the Stokes 
drag (see the lower inset in Fig. 2, where the fit gives 
v =8.4 cS instead of 10 cS). Then by subtraction the 
viscous drag one obtains the reduced wave drag force, 
[F — F v )/F c , as a function of the velocity. This plot 
clearly exhibits the continuous bifurcation to the wave 
drag state at V c = 16.5 cm/sec, that is rather close to 
the theoretical value of 17.1 cm/sec. Here F v , F c are the 
viscous drag force and the critical value of the full drag 
force at the onset, respectively. 

To emphasize a role of a contact line relocation due to 
the ball displacement, two set of data-one with a feed- 
back as in Fig. 2, and another without-are shown in Fig. 3 
for silicone oil of 50 cS. Since the capillary force exceeds 
the wave drag force near the onset, the relocation of the 
contact line due to a ball motion and variation of capil- 
lary forces in this way diminish partially or completely 
the result of the wave drag force in the measurements. 
So, the transition is almost smeared out in the measure- 
ments without the feedback. One should mention that 
in our set-up it was impossible to actuate the feedback 
with a ball less than 1.57 mm dia. Indeed, the viscous 
drag reduces with a radius between linear and square de- 
pendence but the stabilizing magnetic forces reduces in 
a third power. Fig. 4 shows the data of the reduced drag 
force vs the reduced velocity for three different fluids with 
3.14 mm ball: silicone oil(50 cS), glycerol- water (30 cS), 
and glycerol- water (46 cS). It is obvious from the plots 
that the increase in the wave drag force strongly depends 
on fluid properties (y and a). We would like also to point 
out that the wave drag for water(see the inset in Fig. 4) 
increases even more dramatically. The inset presents the 
data for water with 3.14 mm ball in the narrower range 
of the control parameter, V/V c , together with a fit by 
the stationary Ginzburg-Landau equation with a field: 
e£ — a£ 3 + h = o, where £ = (F — F v )/F c is the order 
parameter, e = V/V c — 1 is the control parameter, a is the 
nonlinear coefficient, and h is the field, which plays a role 
of the smearing factor. As we found from the fits both a 
and h depend on v and a of fluids, and ball sizes. The 
scaling found experimentally, is a ~ rja/d 2 . Indeed, using 
this scaling all the data for various fluids and ball sizes in 
the vicinity of the transition collapse onto one curve(see 
Fig. 5 for 5 different sets of the data). Thus, the transition 
to the wave drag state is continuous, and the drag force 
behaves above the threshold as (F — F v ) ~ \JV — V c . 

Fig. 6 presents an image of the water surface waves just 
above the transition that provides information about the 
opening angle of the Cherenkov cone and the wavelength. 
The inset in Fig. 7 shows the dependence of A in water on 
the velocity together with the fit. The main plot presents 
the data for cos 9 as a function of the velocity. Due to 
dispersion the phase velocity of the waves depends on the 
wavenumber. We used the fit of the wavenumbers vs ve- 
locity from the inset in Fig. 7 to obtain the corresponding 
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values of c(k) in eq.(3). The result of the calculations of 
cos 9 according to eq.(3), is shown by a solid line in Fig. 7. 
If one takes into account that no fitting parameters were 
used in this procedure, the agreement with the data is 
rather surprising. Particularly, the critical velocity, V c , 
obtained from extrapolation of 9 — ► 0, is sufficiently close 
to the value, found from the force measurements. 

Thus, the experiments presented shows unambiguously 
that the transition to the wave resistance state is continu- 
ous one, contrary to the theoretical prediction M. What 
can be the reason for this discrepancy? The theory is 
based on the Kelvin model (?]] , which is a reasonable ap- 
proach to long wavelength gravitational waves ^|. In- 
stead of a real ship-like object Kelvin considered a mov- 
ing pressure point applied along its course on a water 
surface [Q. In a case of an object much larger than A c , 
one can neglect a pressure redistribution due to capillary 
effects, and the Kelvin model works extremely well for 
long gravity waves. In the opposite relocation of 

a contact line causes redistribution of the pressure. That 
results in a bump behind the moving object observed 
in experiments. This factor makes the applicability of 
the Kelvin model questionable. The role of capillarity 
and wetting is clearly demonstrated experimentally by 
switching off the feedback control, as we discussed above. 
Moreover, by using a thin wire as an object one can in- 
crease drastically the role of the capillary (wetting) forces 
in a force balance with the viscous and wave drag forces. 
Then strong fluctuations of the contact line position due 
to hysteretic behavior of wetting front, result in strong 
scatter of the data in the vicinity of the transition which 
we actually observed in the experiments with wires of 0.3 
and 0.7 mm. This effect leads the authors of Ref. || to 
an erronous conclusion about the type of the transition 
to the wave drag state. 
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of us (V.S.) is gratefult to M. Assenheimer, H. David- 
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able and helpful advices, particularly on initial stages of 
the research, and to M. Assenheimer and J. Groshaus for 
participation in early versions of the experiment. This 
work was supported by the Minerva Center for Nonlin- 
ear Physics of Complex Systems. 
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FIG. 1. Experimental set-up: FC-channel with a fluid; 
CYL- cylinder to measure a deviation from an initial position 
by eddy current gauge(EMG); IA- amplifier, and PS-power 
apmlifier; C- coils of electromagnets. 

FIG. 2. The drag force vs the velocity for silicone oil 
DC200/10 cS and 3.14 mm ball. The upper inset: the same 
data but in the reduced force. The lower inset: viscous drag 
force below the transition vs velocity. Solid line is the second 
order polynomial fit. 

FIG. 3. The reduced drag force vs the reduced velocity for 
a silicone oil DC200/50 cS and 3.14 mm ball: full squares- with 
feedback, open squares-without feedback. 

FIG. 4. The reduced drag force vs the reduced velocity 
for three fluids and 3.14 mm ball: squares- DC200/50 cS; 
circles-glycerol- water 30 cS; triangles-glycerol- water 46 cS. 
The inset: the same for water, solid curve is the fit by the 
Ginzburg-Landau equation with a field(see text). 

FIG. 5. The scaled data of the reduced drag force vs 
the reduced velocity for 5 different fluids with 3.14 mm 
ball except otherwise mentioned: full squares-water; full 
circles-DC200/50 cS; triangles- glycerol- water 46 cS; open 
squares-glycerol-water 30 cS; open circles-DC200/50 cS, 2.35 
mm ball. 



FIG. 6. Image of the surface waves on water. 
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FIG. 7. cos# of the Cherenkov cone as a function of the 
velocity for the surface waves on water. The solid line is cal- 
culation based on eq.(l) and (3) with use of the data on k(V), 
presented in the inset. 
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